Obtaining cosmic ray propagation parameters from diffuse VHE 
7-ray emission from the Galactic center ridge 



I. Biisching 

North-West University, Potchefstroom Campus, South Africa 

O.C. de Jager 

North-West University, Potchefstroom Campus, South Africa 

J. Snyman 

North-West University, Potchefstroom Campus, South Africa 

ABSTRACT 

The recent discovery of diffuse, VHE 7 radiation from the Galactic center 
ridge by the H.E.S.S. telescope allow for the first time the direct determination 
of parameters of galactic cosmic ray propagation models. In this paper we show 
that the diffuse 7-radiation near the Galactic center may be explained by the 
interaction of VHE cosmic ray (CR) protons with the interstellar gas located 
in several giant molecular clouds leading to a measurement of the cosmic ray 
diffusion coefficient for the galactic center region of k — 1.3 kpc 2 Myr _1 for a 
mean proton energy of ~ 3 TeV, if we assume that the CR protons originated 
from a supernova event (Sgr A East), which took off about lOkyr ago. This value 
of k is ~ 5 to 10 times smaller than the locally measured value. 

Subject headings: diffuse emission, cosmic ray propagation, Galactic center 



1. Introduction 

The High Energy Stereoscopic System of four telescopes offers currently the best angular 
resolution for the study of VHE 7-rays from cosmic sources (Aharonian et al. 2006). With an 
angular resolution of ~ 0.08°, the H.E.S.S. Collaboration was able to resolve 7-rays associated 
with molecular clouds in the galactic center region (Aharonian et al. 2006): Whereas a 
relatively good correlation was found between the 7-ray and CS (with the latter measured 
by Tsuboi et al. (1999)) surface brightness distributions within 150 pc (£ ~ ±1° along 
galactic longitude) from the galactic center, this correlation degraded at a distance of 200 
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pc from the galactic center (i.e. at i ~ 1.5°). A strong indicator that this diffuse component 
is indirectly associated with a source at the galactic center is the similarity of the spectral 
indices of the point source HESS J1745-290 at the center and this newly discovered diffuse 
extended emission. Furthermore, the relatively good correlation within i ~ ±1°, but with 
degrading correlation beyond that, suggests that we are dealing with a source (e.g. the SNR 
Sgr A East or the central Black Hole Sgr A*) at the GC, which was active for some time in 
the past and that we are now seeing the high energetic particles diffusing from this central 
source (Aharonian et al. 2006). The most likely primary species responsible for the gamma- 
ray emission is then protons, since electrons would have to compete against synchrotron 
losses, resulting in a spectral steepening towards large distances. Apart from slight effects 
of energy dependent diffusion, protons, in a good approximation, do not loose energy within 
this environment, resulting in an approximate invariant spectral index with distance. 

Whereas these results are important from an Astrophysical viewpoint, it is also of im- 
portance from a cosmic ray viewpoint, i.e. the study of cosmic ray propagation in our galaxy: 
Aharonian et al. (2006) suggested that the diffusion coefficient for protons in the 4 to 40 TeV 
range should be less than 10 30 cm 2 s _1 (or 3.5 kpc 2 Myr _1 ) as a result of enhanced turbulence 
and higher magnetic field strengths in the GC region. The H.E.S.S. data therefore offer 
a unique possibility to measure the diffusion coefficient in this part of the galaxy and to 
compare with other measurements of the galactic diffusion coefficient. This is in particular 
important as it has been shown (Biisching et al. 2005) that, given SNR are the main sources 
of CR, the widely used method to obtain propagation parameters by fitting secondary to 
primary data is at least tainted, as the CR primary component then shows strong variations 
in space and time. 

In this paper we will model a transient source at the GC with an activity time-scale in 
the past. By solving the transport equation for proton propagation along the galactic plane, 
we will obtain the range of diffusion coefficients, which fit the observed HESS profile best. 
To do this we will also model the gas distribution as traced by the CS emission. 

Tsuboi et al. (1999) was the first to obtain a full coverage of the Galactic Center Bow 
(GCB) and the molecular cloud structures in CS in the region of interest. Since the gamma- 
ray surface brightness is reflected by the line-of-sight integral of the product of the cosmic 
ray and gas densities, we will introduce 3D structures in the GC region (even though not 
strictly modelling the GCB itself), such that the line-of-sight integral through these model 
structures reproduce the true line-of-sight integrals through the gas density within ~ 5% 
accuracy, as described in more detail below. 

Using the diffusive model of CR propagation, together with observations made by the 
High Energy Stereoscopic System (H.E.S.S.) (Aharonian et al. 2006) it is possible to get an 
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estimation of the diffusion coefficient controlling cosmic ray (CR) transport in the galactic 
center region. 



1.1. 7-rays from pion decay 

The omnidirectional (i.e. integrated over solid angle) differential 7-ray source function 
q 7T o(E 1 ,r) at the position f = (l,b,r) for the decay it — > 2j is given by (Busching et al. 
2001) 

qae,^ = 2 r g^fl d7 , (i) 

where 7^ is the pion Lorentz factor. The lower boundary of the integration is given by 

V ~ m n c 2 4 £ 7 ' U 

The pion source function then is 

POO 

Qn»(l^r) = Pgas(r)c / (3a^(j p ,j w o)N p (j p ,r)d-f p (3) 

where a pp is the total cross section for pion production in pp collisions and N p the CR proton 
spectrum. The differential photon flux from the decay of CR induced n°s from the direction 
(/, b) is given by integrating Eq.l along the line of sight 

dN{E 1 ,l,b) 



dt dE y dQ 4tt 



l - J q7r o(E 7 ,r)dr (4) 



The above calculation is for pion production from pp interactions only. The effect of the 
known chemical composition of the ISM can be taken into account by increasing the total 
pion production cross section by a factor of 1.30 (Mannheim & Schlickeiser 1994). 



2. Reproducing the diffusive 7-ray emission from the Galactic center 

For our studies we reproduce the diffusive 7-ray emission calculating the line of sight 
integral Eq. 3 for various (/, b) combinations. As we are only interested in relative intensities, 
we can use instead of Eq. 3 the relative emissivity 

poo 

e(Z,6) oc / p(l,b,r)N CR (l,b,r)dr (5) 

Jr=0 

where p(l, b, r) is the target material density as a function of galactic coordinates and 
Ncr{1, b, r) is the calculated CR density at the point defined by triplet (I, b, r). 
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2.1. 



Gas distribution near the Galactic center 



The inner ±150 pc region of the Galaxy contain interstellar H2 gas of about 2 to 5xl0 7 
solar masses (Tsuboi et al. 1999) in a rather complex setup of molecular clouds. For our 
analysis we assumed that the target material density function can be adequately described 
by the superposition of five spherical Gaussian functions upon an asymmetric Gaussian base. 
These functions represent the molecular clouds associated with the radio arc of Sgr A, Sgr B, 
as well as the longitude varying line-of-sight projection effect of the "Galactic Center Bow" 
as described by Tsuboi et al. (1999). Even though we are able to reproduce the observed 
line-of-sight gas densities within about 5% from the observed values, uncertainties in the 
exact depth distribution along the r-coordinate in Eq. 5, is expected to result in a ~ 50% 
systematic uncertainty in the final estimate of the diffusion coefficient. 



The CR distribution as a function of spatial coordinates is calculated using the diffusive 
model of CR propagation. We calculate the CR density assuming the CR coming from a 
single SNR event 10 kyr ago, accelerating particles for a certain time. As was pointed out by 
Maeda et al. (2002), the SNR Sgr A East has an estimated age of 10 kyr, but also younger 
ages have been given for that SNR (Rockefeller 2005). For our study, we assume that Sgr A 
East has an age of 10 kyr and was accelerating CR on various timescales less than 10 kyr. We 
are interested in the CR distribution near the source and thus neglect the effects of spatial 
inhomogeneities in the interstellar gas density (which are small for CR protons with the 
energies we are dealing with) by assuming a mean gas density, and also boundary effects by 
imposing boundary conditions for the CR propagation problem at infinity. In this case one 
can find a solution for the propagation equation in the literature (Syrovatskii 1959). Using 
the integral in Eq. 5 the emission seen by H.E.S.S. (Aharonian et al. 2006) can be recreated 
for different diffusion coefficients as shown in Fig. 1 (top). As noted by Aharonian et al. 
(2006), the excess counts observed by H.E.S.S. follow the known target material density fairly 
well, except for the region / > 1°. This provides a gauge for approximating the diffusion 
coefficient governing the CR propagation near the galactic center. 

To obtain better counting statistics, we compare the model and observed H.E.S.S. excess 
counts, integrated over Galactic latitude. The excess counts observed by H.E.S.S. then are 
proportional to the integral 



2.2. Cosmic ray distribution near the Galactic center 




(6) 
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where the CR density N C r has been calculated for different diffusion coefficients, as shown in 
Fig. f (bottom). The obtained fit is normalised to the total numbers of excess counts before 
calculating the corresponding x 2 value. The results of the analysis are shown in Fig. 2 below. 
The optimal value for the diffusion coefficient k was found to be 



which is about 40% smaller than the diffusion coefficient estimated by Aharonian et al. 
(2006). We also note that the minimum reduced x 2 is 1-5 (for ~ 25 d.o.f.), which indicates 
that the data are marginally described by the model. Note that we do observe a well-defined 
minimum in Fig. 2, leading to the abovementioned measurement of k. It would be interesting 
to see if more sophisticated 3D modelling of the total gas distribution result in even smaller 
X 2 values. Because of this uncertainty we estimate a systematic uncertainty of ~ 50% on 
the measured value of k. 



In the last section, we have shown that the diffuse 7-ray emission from the Galcatic 
center ridge can be explained by CR hadrons which propagation can be described by a 
diffusion coefficient of 1.3 kpc 2 Myr _1 . This diffusion coefficient is valid for CR generating 
the bulk of the 7-ray emission in the H.E.S.S. energy range. 

To compare this value with that obtained by fitting local CR data, as derived e.g. by 
Moskalenko et al. (2002), Jones et al. (2001) or Maurin et al. (2002), we have to estimate 
the energy of the CR probed here. For this investigation, we approximate the sensitivity 
of the H.E.S.S. telescope array by a box function from 0.4 TeV to 20TeV. For the CR 
proton spectrum we assumed a power law with the observed photon index Y = 2.29 ± 0.27 
(Aharonian et al. 2006). We used the Pythia (Sjostrand et al. 200fa,b) event generator 
package to calculate the number of 7-rays with £" 7 > 0.4 TeV, for different CR proton 
energies. The result of this calculation is shown in Fig. 3. The maximum number of photons 
is produced by CR protons with energies of about 2.2 TeV, given a proton spectral index of 
2.29. Note that the location of the maximum depends on the proton spectral index. Given 
the uncertainties stated by Aharonian et al. (2006), we find that CR protons with energies 
in the range 1.7 TeV to 3 TeV contribute mostly to the 7-rays seen by H.E.S.S. Assuming a 
diffusion coefficient of the form 



where £ = 1 GV and £ the particle rigidity we find k = 0.013 kpc 2 Myr 2 , which is signif- 
icantly smaller than the values found by fitting local CR data. However, k for the latter 



k = f.3kpc 2 Myr 



-1 



(7) 



2.3. Mean CR energy 




(8) 
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Fig. 1. — Calculated emission sky maps (top) together with the calculated excess counts (bot- 
tom) shown as the solid lines. The histogram indicate the excess counts from the H.E.S.S. 
observation (Aharonian et al. 2006). The figures (a), (b) and (c) where generated for diffu- 
sion coefficients with values 0.3 kpc 2 Myr~ l , 1.3 kp(?Myr~ x (best fit) and 15.1 kp(?Myr~ x 
respectively. 
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Fig. 2. — Reduced x 2 values plotted for different diffusion coefficients. We find a well defined 
minimum for a diffusion coefficient of k = 1.3 kpc 2 Myr~ l . 
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Fig. 3. — Number of photons with E 1 > 400 GeV for different CR proton energies assuming 
a CR proton spectral index of s = 2.29. 
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range from 0.0535 kpc 2 Myr" 2 to 0.201 kpc 2 Myr" 2 , as compiled by Maurin et al. (2002) (with 
references therein) . This finding can be well explained by enhanced turbulence and a higher 
field strength of the interstellar magnetic field in the Galactic center region. We note how- 
ever that the rigidity dependence of 0.6 can only apply to a limited energy (rigidity) range, 
since k must always be larger than the Bohm limit. 

3. Summary and Discussion 

We have shown that the progress in the imaging Cherenkov technique now makes it 
possible for the first time to measure the CR diffusion coefficient in other parts of the 
Galaxy. A diffusion coefficient of k, ~ 1.3 kpc 2 Myr _1 appears to be well measured from the 
data, although we have to add a ~ 50% systematic undertainty arising from uncertainties 
in the actual 3D gas density distribution, as well as uncertainties on the epoch when the 
central source activity started, which was assumed to be 10 kyr in this paper. The most 
likely central source is Sgr A East, for which the epoch of onset of activity is ~ 10 kyr, but if 
the central source was the central massive black hole Sgr A*, the timescale of activity would 
be much less certain, resulting in even larger uncertainties on k. Note however that, given 
an initial epoch of onset of activity (i.e. 10 kyr), estimates of k appear to be fortunately 
robust against uncertainties in the detailed time profile of particle acceleration within 5 kyr 
after the SN explosion. For example, reducing the "on" -time of the source by a factor of 
five, decreases the diffusion coefficient by 40%. Such detailed studies will however be treated 
in a subsequent paper. Thus, if Sgr A East was the source of CR, then we have a relatively 
accurate masurement of k. 

In general, following the notion of Aharonian et al. (2006) to constrain the cosmic ray 
diffusion coefficient k from the spatial distribution of diffuse 7-rays resulting from impulsive 
injection of a CR source at some time in the past, we have shown that one can obtain unique 
measurements of the diffusion coefficient, provided that the gas density distribution, as well 
as the epoch of onset of central source activity is known. We also showed that the CR 
diffusion coefficient in the Galactic center region is significantly smaller than that obtained 
by fitting local CR data. Our understanding of turbulence theory is however still too limited 
to understand how diffusion coefficients scale with the turbulence SB, the correlation length 
associated with this turbulence, the total magnetic field strength B and rigidity dependence 
for perpendicular and parallel diffusion. Hopefully our new measurement of k will help to 
constrain results from turbulence theory. 
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